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1. Introduction. 

When we quantise a field in curved space, many interesting phenomena arise. For 
example in an expanding Universe we have particle creation, which means that a spacetime 
which looks empty in terms of particle modes natural in the past, may look full of particles 
on using co-ordinates natural to the future. The Minkowski vacuum appears to have a 
particle flux for an accelerated observer. For black holes, imposing vacuum conditions 
at past null infinity gives Hawking radiation in the future, due to the time-dependent 
gravitational field of a collapsing object. and references therein.) 

For a complete description of the physics, we need also the backreaction of the quan- 
tum field on the geometry. At this point we need a consistent theory of quantum gravity, to 
deal with the gravity loops which arise. For this reason backreaction questions have been 
traditionally restricted to semiclassical order (apart from recent studies on low dimensional 
gravity theories [H). 

Strings provide a consistent theory of quantised matter and gravity, so we would like 
to investigate the issue of particle creation and backreaction in string theory. Several 
examples of strings in time-dependent backgrounds were studied by de Vega et. al 0. 
In their first quantised treatment of the string, they find that an existing string may 
get excited to higher mode levels in moving through the time-dependent geometry (e.g. 
spacetime with a gravitational wave), but it remains one string. Thus at their tree level 
calculation we do not see the creation of a 'flux of strings', which would be the analogue 
of the 'creation of particles' for field theory in curved spacetime. 

In what way do we expect to see 'created particles' in a first quantised language? Sup- 
pose a Universe starts in the vacuum, and undergoes expansion which creates particles. 
Consider the two point propagator computed for spacetime points in the region where the 
'bath' of created particles is nonvanishing.We may compare this situation to the computa- 
tion of the propagator at finite temperature, where also a bath of ambient particles exists. 
In the latter case we know that the propagator depends on the temperature, and with this 
modified propagator Feynman diagrams can be computed to develop perturbation the- 
ory. With our 'created particles' there should be a similar modification of the propagator, 
except that the 'bath' would be neither exactly thermal nor time-independent in general. 

In this paper we take the simple case of a scalar field, and show how a careful han- 
dling of the 'sum over paths' in the first quantised propagator allows us to obtain the 
modifications needed for incorporating the effects of created particles. If we compute the 
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analogue of the Polyakov path integral for a particle, we get a two point scattering ampli- 
tude between the 'in' and 'out' vacuum states. Such a naively computed first quantised 
propagator does not directly 'see' the created particles. What we need instead is a prop- 
agator that reproduces the causal amplitude obtained by taking the expectation value of 
field operators in a specified state. Our 'first quantised' calculation reproduces this propa- 
gator. A density matrix formalism develops naturally, and the particle fluxes produced in 
time dependent geometries become an essential part of the first quantised description. It 
appears that only such an extension of the Polyakov path integral can adequtely describe 
the physics of quantised gravity plus matter. 

For the string case we have the further issue of obtaining the background field equa- 
tions by a /3-function calcaulation. The tree level /9-function gives the classical vacuum; 
but we need the higher string loop contributions to take into account the backreaction 
of created particles. We argure that our modified propagator, rather than the Polyakov 
prescription, should be used in computing the loop amplitudes. 

The issue. 

Consider the free massive scalar field. For the 'propagator' from one spacetime point 
a;i to another spacetime point X2 , the analogue of the Polyakov path integral gives 0] 



where r parametrises the world line, and ^(0) = xi, X{1) = X2- The action is multiplied 
by i because we are in spacetime with Minkowski signature. The regulator term — eA 
ensures convergence of the path integral. 

What precisely is the answer we obtain? In spacetimes with Euclidean signature there 
is a unique Green's function, determined by the boundary condition that the amplitude 
to propagate to infinity falls to zero in all directions. Such is not the case for spacetimes 
with Minkowski signature, where various Green's functions can be defined, differing from 
each other by solutions of the homogeneous field equation. For describing the choice of 
boundary conditions implicit in ( |1.1| ) let us assume that the spacetime does not expand 'too 
rapidly' in the far past {t — oo) or in the far future {t —>■ oo). Then we obtain naturally 
defined vacuum states for the Klein-Gordon field, |0 >in and |0 >out for t —oo and 
t — i> oo respectively. For a time-dependent geometry, in general |0 >out7^ |0 >in- In this 
situation it is known that the RHS of ( p..l| ) computes 
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which is the analogue of the Feynman propagator in the time-dependent geometry. 

Such 'in-out' amphtudes are relevant for computing scattering amplitudes where a 
finite number of 'in' particles scatter to a finite number of 'out' particles. By contrast, 
for an 'initial value problem', with the initial state the 'in' vacuum, we would compute 
amplitudes like in < 0\T[(j){z2)<p{zi)]\0 >in- In particular, we know that for a source term 
to Einstein's equations we should use 'true expectation values' < 'i/'|T^;y|'i/' > for the stress 
tensor, not 'in-out' vacuum amplitudes |^. We illustrate (in section 5 below) the difference 
between these two different kinds of amplitudes by a simple example in 1 -|- 1 spacetime. 
We compute the spatial average of in < Q\T^jj{r], x)|0 >in (called < T^^{r]) >in in) and of 
out < 0|T^,,(?7,x)|0 >in /out < 0|0 >in (called < T^^{r]) >in out) in a Universe which is in 
the initial vacuum and undergoes sudden expansion at = 0. (77 is the conformal time.) 
For 7] < 0, < Too >m in IS just the vacuum energy. For 77 > it becomes the vacuum 
energy plus the energy of created particles. But < Too >in out is just the vacuum energy 
in both ranges of ry. Worse, < Tn >in out is complex, so it cannot be a source term for the 
gravitational field. 

To compare 'in-in' and 'in-out' amplitudes we note that in < Q\ = Cq out < 
0| + C2 out < + C4 out < 4| -|-. . . where < 2| is a state with one pair of 'out' particles 
etc. The first quantised path integral described above computed only 'in-out' amplitudes. 
This may not appear to be a serious shortcoming of the first quantised approach in the 
particle case because one could compute successively out < 0|0 >im out < 2|0 >im etc. 
and thus reconstruct the result of the initial value problem. But this is not an adequte 
treatment for the string case, which was the motivation for studying the first quantised 
approach in the first place. With strings we do not start with a Lagrangian governing the 
field; instead we discover backgrounds satisfying the field equations by demanding that 
string propagation in the background be 'consistent'. But which propagator should we 
demand consistency for: the 'in-in' one or the 'in-out' one? These two choices would in 
general give different answers for the background field (considering f3- functions to one loop 
or beyond). The discussion above indicates that the 'in-in' propagator is the one which 
correctly gives the stress tensor of the created fiux of 'out' particles, and is therefore the 
one that should be involved in the multiloop contributions to the /3-function. 

Can we develop a natural extension of the Polyakov prescription so that we compute 
'in-in' amplitudes? In the second quantised field language, 'in-in' perturbation theory is 
described in the 'real time formalism' of many-body theory. We outline this formalism 
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below, and then give our results on how this formalism arises naturally and covariantly in 
the first-quantised 'sum over paths' language. 

Note: In this paper we use the term 'Polyakov path integral' for any prescription where 
one sums over spacetime co-ordinates on the world sheet /world line and the metrics on this 
world sheet/world line. This approach was developed primarily for Euclidean target space 
metrics. We are assuming here that the prescription of summing over co-ordinates and 
metrics has been extended in some way to spacetime with Minkowski signature. Thus the 
world sheet will also have Minkowski signature, and construction of higher genus surfaces 
may require explicit use of interaction vertices. 

The real time formalism. 

To obtain a perturbation theory for the initial value problem, one uses the 'real time 
formalism', developed in the context of time-dependent many body theory Suppose 
we start in the vacuum |0 >in at t = — oo. We wish to compute the expectation value of a 
product of observables Oi, defined at times ti, in the state |0 >in- This expectation value 
can be expressed as TrjpoOi . . . On}/Trpo, with po = |0 >in m < 0| the density matrix for 
a pure state. Because both the bra and the ket states composing po are given at t = — oo, 
the interaction Hamiltonian acts along the 'time path' leading from t = — ootot = oo, 
and then back to t = — oo. Correspondingly, we need propagators from any time on the 
first or second time path segment to any other time on either path segment. We collect 
these four propagators (appropriate to the four possible choices of time path segments) 
into a 2x2 matrix propagator. Vertices arising from the interaction Hamiltonian acting 
on the first time path segment have the usual coupling —ig, while those on the second 
path segment have coupling ig. With this extension of rules, the computation of Feynman 
diagrams gives correlation functions for the density matrix po, inserted at t = — oo. 

Note that the above formalism computed 'in-in' correlation functions, regardless of 
whether or not the 'in' vacuum differed from the 'out' vacuum. The density matrix po can 
be replaced by an arbitrary density matrix p. A special class of p emerge in the development 
of perturbation theory: density matrices expressible as 'exponential of quadratic in the 
field operator'. For such p the Wick decomposition holds: the expectation of a product of 
several operators decomposes into sums of products of expectations for pairs of operators. 
The propagator of the real time formalism encodes a choice of such a 'special' density 
matrix. Deviations of p from this special class give 'correlation kernels', which are handled 
perturbatively just like interaction vertices. 
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The relevance of the real time formalism for field theory in time dependent geometry 
is easy to see. Firstly, cosmology appears set up as an initial value problem, rather than 
a scattering problem. Secondly, 'particle creation' is generally obtained as a Bogoliubov 
tranformation on a vacuum state. The particle flux produced this way is described by a 
density matrix of the form 'exponential of quadratic in field'. This flux is therefore pre- 
cisely of the form that can be incorporated into the propagator of the real time formalism, 
so that backreaction computations can be examined by loop calculations using the matrix 
propagator. Thirdly, the density matrix of the real time formalism can be taken to de- 
scribe a mixed state, whereupon one can study correlations in an evolving distribution of 
particles. Note that the 'periodic imaginary time' formalism for finite temperature corre- 
lations applies only to time independent thermal distribution functions. But in a theory 
with gravity, a nonempty distribution of particles gives rise to a nonvanising stress tensor, 
which leads to a time dependence of the spacetime geometry in general, and a correspond- 
ing evolution of the particle distribution. Thus the real time formalism must be used in 
developing kinetic theory for a theory with gravity; the imaginary time formlism can only 
provide quasi-static approximations. 

The above description of the real time method pertained to the second quantised field 
language. In order to apply the method to strings, we need to obtain it in a first quantised 
path integral language. At first one might think that this would not be possible, since a 
first quantised path integral seems to describe the propagation of just one particle, while 
the real time formalism is built to handle fiuxes. Nevertheless, it turns out that one does 
find the entire structure of the formalism, when one carefully quantises the action of, say, 
a scalar particle. The phrase 'first quantised' is not really correct for the resulting path 
integral; what one gets is a 'proper time representation' for the real time propagator.0 
The fact that one reaches this object after starting with what looks like just one particle, 
refiects the ambiguity of the particle concept in curved space. The Polyakov path integral 
(and its particle analogue) chooses a prescription to define the propagation amplitude such 
that this ambiguity is supressed. In restoring the ambiguity, one regains the freedom to 
have an 'in-in' formulation, and also the choice of an initial density matrix. 

Approach and results. 

We start with the geometric action m J ds for a scalar particle moving in a spacetime 
with Minkowski signature, ds = (dX^dX^)^/"^ has a sign ambiguity. Along a trajectory 

^ The author thanks L. Ford for pointing out this issue of terminology. 
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contributing to the sum over paths, the path can change several times between being 
timehke and being spacehke. At each change dX^^dX^ passes through zero, and we have 
to make the choice of root afresh. We might try to evade the problem by passing to the 
quadratic form of the action. One often thinks of this form of the path integral as being 
a sum over functions X^{t) and over metrics grrij) on the world line of the particle. In 
fact this would be the analogue for the scalar particle of the Polyakov prescription. But in 
constraining the square root action to bring it to the quadratic form, what we really get 
is a sum over the X^{t) and over the einbein A(r) on the world line. Since the einbein is 
the square root of the metric, two signs of the einbein correspond to the same (7tt(t); this 
appears to be a reflection of the sign ambiguity in the original square root action. 

To evaluate the path integral wc gauge flx A(r) to a constant by using the freedom 
of diffeomorphisms along the world line. If the world line geometry were described by the 
metric grrij), it is evident that the only degree of freedom remaining of this geometry 
would be the length A of the entire world line, and A should be integrated from to oo. 
But since we have instead a sum over the einbein A(t), the situation is more complicated. 
Diffeomorphisms can fix A to a positive constant in any interval where A is positive, and 
to a negative constant in any interval where it is negative. These two kinds of intervals 
will alternate along the world line for a typical realisation of A(t) in the path integral. 
We can perform the integral over the value of the constants Aj obtained in these intervals, 
with the integral running over positive or negative reals respectively for the two different 
kinds of intervals. To complete the evaluation of the path integral, we need to know what 
amplitude to assign to a crossover of A from positive values to negative or the other way 
round. 

This amplitude is not supplied by the action, and so is supplemental information 
needed to make sense of the propagation amplitude. It turns out that the freedom of choice 
here is exactly that of a density matrix of the 'special' kind mentioned above (exponential 
of quadratic in the field) and the choice of a 'time path'. These are precisely the ingredients 
which define the matrix propagator of the real time approach. The two different signs of 
the einbein A on the world line correspond to the two segments of the time path of the real 
time formalism, and so furnish the two dimensional vector space on which the 2x2 matrix 
propagator acts. 

More specifically, we can describe the first quantised propagator in the following terms. 
We consider both the propagation e~*('-'+™^) and the propagation 6*^'-'+'"^) on the world 
line. We may collect these two possibilities into a 2x2 matrix form, getting a world line 
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Hamiltonian diag[(D + m^), — (□ + m^)]. Near the mass shell □ + = the path integral 
needs regulation. If we add —iti to the Hamiltonian (/ is the 2x2 identity matrix) then 
we get a diagonal matrix propagator, with ( |1.2| ) and its complex conjugate as the first 
and second diagonal entries. But if we regulate instead by adding — ieM, where M is a 
non-diagonal matrix, then we get the matrix propagator for perturbative calculation of 
Green's functions in an 'exponential of quadratic' particle flux. 

We compute M for a thermal distribution in Minkowski space; for the 'closed time 
path' mentioned above, and for the time path used by Niemi and Semenoff [^. These 
two M matrices are not the same, which reflects the fact that the matrix propagator 
depends on both p and the choice of time path. For a curved space example, we consider a 
1 + 1 spacetime with 'sudden' expansion. We compute M to obtain the matrix propagator 
appropriate to the initial value problem with state the 'in' vacuum; i.e. for p = in|0 >< Q\in 
and time path beginning and ending at past infinity. 

We conclude with a discussion of the significance of our results for strings, which were 
the motivation for this study of the first quantised formalism. The usual first quantised 
path integrals for strings would give the analogue of (|1.2|) , which corresponds to specific 
boundary conditions at spacetime infinity. To handle phenomena involving particle fiuxes, 
we propose extending the world sheet Hamiltonian to a 2x2 Hamiltonian as in the above 
particle case. Thus we would allow not only classical deformations of the background 
field (analogous to the 'exponential of linear' p in the particle case) but also particle fiux 
backgrounds (corresponding to 'exponential of quadratic' density matrices). Studying (3- 
function equations [§] for this extended theory should give not equations between classical 
fields but equations relating fields and fiuxes. The latter kind of equation, we believe, 
would be natural to a description of quantised matter plus gravity. 

The plan of this paper is as follows. Section 2 reviews finite temperature perturbation 
theory, and discusses the significance of 'exponential of quadratic' density matrices for the 
curved space theory. Section 3 translates the finite temperature results of fiat space to first 
quantised language. Section 4 gives a curved space example. Section 5 discusses strings. 
Section 6 is a summary and discussion. 
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2. Propagators in the presence of a particle flux. 

2.1. Review of the real time formalism. 

The real time method has been extensively developed in recent years, in a variety of 



contexts. A good description of the closed time path and its meaning is given in [|T0[, jTl 
We give a brief summary below. 

Consider a scalar field in Minkowski spacetime (metric signature -\ — ■ ■ ■ — ) 



S = j dx{\d^^d>^<j) - \m'c\? - ^<^4) (2.1) 



We have the following relation for Green's functions [12 



< 0|T{(/>(a:i) . . .(/>(a:O}|0 >=< 0|t/-\oo)T{,^o(xi) . . . c\)^{x^)e ' ^*'^}C/(-oo)|0 > 

(2.2) 



where 



U{t) = Te"'-^*'=— (2.3) 
(P{x) = U-\t)Mx)Uit) (2.4) 



and (f)Q is the free field. 
From ( |2.3| ) we have 



C/(-cx))|0 > = |0 > (2.5) 

If the vacuum is assumed stable: 

< 0\U-\oo) = e'^ < 0| (2.6) 
then we obtain the usual perturbation series 

^ni™ ^ nin^ < O|r{c/>o(xi) . . .M^n)e-' ^^--''''''-'^''^}\0 > 

< e Jt'=-oo U > 



If the vacuum is not stable (i.e., ( p.6| ) does not hold), then we get interaction vertices 
from U~^{oo) in ( pj.2| ). The interaction Hamiltonian acts from t = — oo to t = oo because 
of the exponential in ( |2.2| ) and from t = oo to t = — oo because of U~^{oo). The Green's 
function can be written as 



< 0|T{<^(xi) . . .(^(x^llO > = < 0\Tc{Mxi) ■ • •<^o(^n)e ^ /c ^*'^-(*')}|0 



> (2.^ 
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where H[^t is integrated along the time path C running from t = — ootot = oo and back 
to t = — oo. Tc represents path ordering along the path C. 

We label field operators on the first part of C with the subscript 1 on the second part 
of C with the subscript 2. Thus for the two point function of the scalar field there are four 
possible combinations of subscripts, and the corresponding propagators are collected into 
a matrix: 

^ \ - ( 0|T[(/)(X2)(/'(X1)]|0 < 0J(/)(X1)(/)(X2)|0 \ , . 

^ V .n<0|(/>(x2)(/>(xi)|0>,„ ,„<0|T[</>(x2)(/>(xi)]|0>,J ^^-^^ 

where we have added the subscript 'in' to the vacuum to emphasize that it is the vacuum 
state at t = — oo. (T describes anti time ordering.) 

Perturbation in the coupling A is described by the following diagram rules [0. The 
external insertions are all of type 1. (This can be seen from the way we arrived at ( |2.8| ) 
from (|2.2|) .) The vertices from the perturbation term must have either all legs of type 1 
or all legs of type 2. Attach a factor {—iX) to the vertex in the former case, [iX) in the 
latter. The propagator is given by the matrix (p.9| ) and can connect vertices of type 1 to 
vertices of type 2. Summing all Feynman diagrams with these rules gives the correlation 
function for the initial value problem with initial state the vacuum at t = — oo. 

As a simple example where the second leg of the time path makes a difference, consider 
the free scalar field in + 1 spacetime dimensions. Let the only perturbation be the time 
dependent term ^ f dtd{t — to) '■ (f^if) '■■ We compute to first order in ^ the two point 
correlator for field insertions at ti < to < ^2, using time path C: 

< (I){t2)(j){ti) > = 2fi[< 0\(P{t2)(Pito)\0 >< 0|(/)(to)<^(ti)|0 > 

(2.10) 

+ < o|(^(to)0(t2)|o X omo)(p{h)\o >] 

Without the second leg of the time path, we just get the first term on the RHS. The 
latter answer gives the amplitude to start in the interaction picture vacuum at t = — oo 
and to also end in interaction picture vacuum at t = oo. The missing term corresponds 
to correcting the state at t = oo away from the vacuum |0 >, to be such that it evolves 
from |0 > after suffering the perturbation at to. Thus the complete answer ( p. 10 ) gives a 
Heisenberg picture expectation value, and is thus the only physically meaningful object to 
be computed in this situation. 
We can rewrite ( p.8|) as 

<0|T{(^(xi)...(^(x^)}|0> = TV{poTc{(/>o(xi)...(^o(xOe-*/c^*'^'-(*')}}/Trpo (2.11) 
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with p = |0 >< 0| the density matrix of a pure state. But we can immediately extend 
the above relation to to the case where po is replaced by an arbitrary density matrix p, 
not necessarily representing a pure state, and not necessarily describing a distribution in 
equilibrium. To develop diagrammatic perturbation theory, however, we need to have the 
Wick decomposition for expectation of operator products. This is achieved for a special 
class of density matrices, which have the form 'exponential of quadratic in the field operator 
Departures of p from this special form are handled perturbatively just like interaction 
vertices. We discuss the significance of this special class of p in the next subsection. 

As an example of mixed state density matrices, let p at t = — oo be of thermal form 
with temperature Then the matrix propagator is 



^^^^ V 27r[n(p)+%oMp2-m2) + 27m(p)5(p2 _ ^ I - J 

where 

n{p) = (e^lPol - (2.13) 

is the number density of particles for the free scalar field at temperature {I3)~^ . With this 
change of matrix propagator, perturbation theory has the same rules as for the vacuum 
case. 

Another approach leading to a similar diagrammatic structure was given by Niemi 
and Semenoff In the periodic imaginary time description of temperature, we evolve 
the field theory in imaginary time from t = to t = — and identify these two time 
slices. In the 'real time' approach of Niemi and Semenoff we use instead the following path 
Ci in complex t space to connect these two slices 

Ci : / : — oo oo 

// : oo^oo-z/3/2 

(2.14) 

/// : oo - i[3/2 ^ -oo - z/3/2 

IV : -oo - i(3/2 ^-oo-i(3 

The concept of time ordered correlation functions is replaced by 'path ordered' correlation 
functions, with t running along the above path Ci. One argues that the parts //, IV of 
Ci can be ignored. Field operators on part / are labelled with the subscript 1 while on 
part /// are labelled with the subscript 2. The matrix propagator becomes 

/ +27rn(p)^(p^-m^) 2n[{nip)inip) + l)]^/^Sip^ ~ m^)\ 

^^^^^ V27r[n(p)(n(p) + l)]i/25(p2_^2) + 27vn{p)d{p^ - m^) J 

(2.15) 
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Perturbation theory based on the contours C and Ci are not equivalent, even for zero 
temperature [0. The propagators ( p.l2| ) and ( p.l5|) are not the same, for /3 — > oo. Thus 
the natrix propagator depends both on the particle distribution and on the time path. 

Note that the temperature dependent terms in ( |2.12| ) or ( |2.15|) are all on shell. This is 
a manifestation of the fact that these terms describe a flux of real particles. The expression 
5{p^ — m?) does not have a Euclidean counterpart. (In Euclidean space one would consider 
{p^ + m?) which is positive definite.) Thus temperature is a phenomenon of Minkowski 
signature spacetime, where Green's functions are subject to the addition of solutions of 
the homogeneous field equation. 

To understand the origin of the on shell terms, let us consider the element Dn in 
( p.l5|) . We can decompose the free scalar field in Minkowski space into Fourier modes. 
Each quantised mode is a harmonic oscillator with frequency c<;(p) = (p^+m^)^/^. The two 
point function < T[^(x2)(^(xi)] > reduces, for each Fourier mode, to < T[q{t2)q{ti)\ > for 
each harmonic oscillator. Thus we focus on a single oscillator (we suppress its momentum 
label). Let t2 > ti. Then 

< qit2)qiti) > = ^e-^("+i)'" < n\T[q{t2)qiti)]\n > / ^ e-^('^+^)'" 

(2.16) 

^ <n+l > < n > ^i^(t2-u) 

The first term on the RHS corresponds to the stimulated emission of a quantum at ti with 
absorption at t2- The second term corresponds to the annihilation at time ti of one of the 
existing quanta in the thermal bath, and the subsequent transport of a hole from ti to 12, 
where another particle is emitted to replace the one absorbed from the bath. The Fourier 
transform in time of (|2.16|) is 

/ dte-'^' < T[q{t2)qiti)] > = -75 5 + Stt < n > 6{u;'^ - u;'') (2.17) 

which gives the matrix element Dn in ( |2.15|) . 

Thus the correction to < T[(^(x2)(^(a:i)] >= Dii{x2iXi) (and other elements of D) 
due to the particle flux is not an effect of interactions with the particles of the bath. This 
correction arises from the possible exchange of the propagating particle with identical real 
particles in the ambient flux. Using such a corrected propagator with the interaction 
vertices gives the interaction of the bath particles with the propagating particle. With this 
understanding of the propagator we see that the real-time formalism is really a covariant 
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concept. In a general spacetime, single particle wavefunctions are just solutions of the 
free field equation. If there is more than one particle populating the same wavefunction, 
then we must keep track of 'exchange terms' in computing correlation functions. The 
matrix propagator does this bookkeeping; the terms with 5{p'^ — m^) in the above would 
be replaced with the appropriate solutions of the wave equation in the spacetime. 

2.2. 'Exponential of quadratic' density matrices. 

When doing perturbation theory on flat space, we usually assume that the vacuum is 
stable: |0 >out= |0 >in (upto a phase). As we saw above, when we do not have this stability 
then we have to use a 2x2 matrix propagator. Stability may be lost either because the 
spacetime geometry is time-dependent, or because there exists a distribution of particles 
in the initial state. With gravity we would typically encounter both these causes: the time 
dependent geometry creates particle fluxes which evolve on (and modify) the geometry. 
Thus it becomes natural to construct the theory with arbitrary initial state and to not 
consider a 'vacuum theory' at all. 

In perturbation theory we deal easily with a special class of states: the coherent states. 
These have the form e^t''^l|0 > where A is an operator linear in the fleld variable (f). The 
operator A generates a change in the classical value of the fleld; shifting the fleld variable 
by this classical value removes the effect of A from Feynman diagrams. In studying flrst 
quantised propagation about a background, this shift has already been made (to zeroth 
order in the coupling). 

The discussion of the real time method revealed that there is another kind of fleld 
conflguration that is also easy to deal with. This is the class of conflgurations expressed 
as density matrices p = e^^'^\ where B is quadratic in the fleld operator. The special role 
of these 'exponential of quadratic' density matrices is due to the fact that Wick's theorem 
extends to correlators computed with such p [^. Thus for operators Ai linear in the fleld, 

< Ai...Ar,>= ^TT{pAi...Ar,} = < > ... < Ai„_,A,„ > (2.18) 

permutations 



We sketch a proof of (|2.18|) in the appendix. 

The physical importance of 'exponential of quadratic' density matrices is seen from 
the different instances where they occur: 

a) The thermal distribution has such a density matrix: 

Pthermal = \{ C"^'^^"!"^ (2.19) 
modes i 
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b) In time dependent geometries one considers the Bogoliubov transformation re- 
lating different expansions of the field variable (we supress the mode index): 



O'infin + 0, in fin 



t 

(2.20) 

out 



(p — a out font + a\utfi 

The 'in' and 'out' vacuua are given by 

Om|0 >in = 0, aout\^ >out = 0. (2-21) 

The 'in' vacuum can be expressed as 

|0 >^n = e'^^-^'-HO >out (2.22) 

so that it contains a flux of 'out' particles. If the initial state is |0 >in then the density 
matrix for the initial value problem is po = |0 >i„ in < 0|. po can be expressed as 

(2.23) 

=lim^^oo [p/Trp] 

where 

p = e''"°"'"°"'e~^"°''*"°"*e^*"°"'"°"' (2.24) 

and the limit is established by inserting complete set of 'out' number operator eigenstates 
between the exponentials in (|2.23| ). 

Thus the 'exponential of quadratic' density matrices 

p = lie' e »«e^' (2. 25) 

i 

unify the treatment of thermal fluxes and the fluxes created by gravitational flelds. An 
ensemble like ( |2.iy| ) might describe the distribution at the initial state of the Universe. In 
fact the notion of 'thermal' distributions has to be extended to the class ( |2.25|) to be useful 
in curved space. Suppose we choose a time co-ordinate t and start with a distribution 
e~^^, thermal for the Hamiltonian giving evolution in t. As the Universe expands, the 
distribution will not remain thermal, in general. Redshifting of wavelengths gives an 
obvious departure from thermal form if the fleld has a mass or is not conformally coupled. 
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But even a massless conformally coupled field departs from thermal form if the time co- 
ordinate t is not appropriately chosen [^. However, the density matrix remains within 
the class (|2.25| ), if it starts in this class, even for a massive field. 

The Bogoliubov transformation is important in obtaining the stress energy of particles 
created in the gravitational field. For example in the black hole geometry the state is the 
Kruskal vacuum. The Bogoliubov transform to the Schwarzschild modes allows us to see 
readily the < T^i, > of emitted radiation at spatial infinity. The class (|2.25|) is closed under 
Bogoliubov transformations relating different field expansions like ( |2.2C1| ). It is the natural 
class to arise in a covariant formulation where the density matrix po discribed above is one 
of the allowed initial state specifications. 

How do we compute the backreaction on the gravitational field from particles described 
by (|2.25| )? For a coherent state, the stress tensor is given, ignoring quantum corrections, by 
the classical value of < T^^, > computed for the classical field value implied in the coherent 
state. But the density matrices ( p.25| ) do not approximate a classical configuration; indeed 



the dispersion in T^^, in each mode is always of the order of T^j^ itself |T^. Can we take 
into account the backreaction from such distributions? By incorporating the 'exponential 
of quadratic' density matrices in the propagator the real time formalism allows us to do 
precisely that. Thus even though such configurations are far from classical the Feynman 
diagram technique can handle them, provided we make the extension to the 2x2 matrix 
propagator. 

What if p was not well approximated by the form e^^^^^-^^^^l No easy route appears 
available for studying backreaction in that case. Small deviations from the above form can 
be handled perturbatively as 'correlation kernels', which behave as extra, nonlocal vertices 
for the perturbation series. For a description of the expansion in correlation kernels and 



interaction vertices, see 14 



Thus in the real time method we are relying on the assumption that for the situations 
of interest the state is well approximated by coherent and 'exponential of quadratic' parts. 
From the discussion above about particle creation, we conclude that allowing coherent 
states alone will not be adequate for a theory of gravity, particularly in regions of strongly 
varying garvitational field. But if we set all couplings (other than gravity) to be small, 
we should be able to validate the above assumption about p, and carry out real time 
perturbation theory. 

We should distinguish two different limits in which the physics of fluxes may be stud- 
ied. One limit is where the collisions are so rapid that approximate thermal equilibrium 
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is maintained at all times, and we need only let /3 be a function of time. The other limit 
is that of kinetic theory, where we assume that collisions are rare; particle wavefunctions 
evolve on the time-dependent background, and collisions between these particles are taken 
into account by perturbation theory. Our approach assumes the latter limit. 



3. The first quantised formalism. 

3.1. The regulator matrix. 

The Feynman propagator for a scalar field in Minkowski space can be written as 

Df{p) = = / dAe^^(p'-"^'+^^) (3.1) 



( pT|) can be used to express Gf{p) in a first quantised language, with p = — U. (See for 
example |]1^, [0.) The Hamiltonian on the world line is □ + m?, evolution takes place in 



a fictitious time for a duration A, and this length A of the world line is summed over all 
values from to oo. 

Let us write the matrix propagator ( p.l5| ) in a similar fashion 

/>oo 

Dns{p) = / rfAe-^^^-^^^ (3.2) 
Jo 

where 



(p2_^2) X l + 2n{p) -2^n{p){n{p) + l) 



{p^-m^)J' \-2y/n{p){n{p) + 1) 1 + 2n{p) 

(3.3) 

{n{p) is given by ( ^.13) ).) This matrix world line Hamiltonian has the following structure. 
If we forget the term multiplying e, then the a = 1 component of the vector state on the 
world line evolves as e*^*-^ ""^ while the a = 2 component evolves as e"*^*-^ \ To 
define the path integral we need the regulation from e at the mass shell p'^ — m? = 0. 
But the regulator matrix is not diagonal, for nonzero temperature. Thus transitions are 
allowed from state 1 to state 2. As we shall show below, in the limit e ^ 0"*", which we 
must finally take, these transitions occur only on the mass shell. (The absolute values of 
the entries in M are not significant, because e goes to 0, but the relative values are.) 

Similarily we can write the matrix propagator ( p.l2|) in the form ( p.2|) with H as in 
but 



^ 1 + 2n{p) -2^n{p){n{p) + l)e-fP^/^- 

' -2v/n(p)(n(p) + l)e^Po/2 1 + 2n{p) ' ^ ' 
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Again we have the a = 1, 2 states propagating on the world hne with Hamiltonians =f(p^ ~ 
m^). Transitions between a = 1, 2 are again of order e, but are different from those in 

Looking at these examples the following picture emerges. To obtain the matrix propa- 
gator in a many-body situation (with 'exponential of quadratic' density matrices) we need 
to consider both evolutions e~'^^^ and e*^^ on the world line. A damping factor is needed 
to define the first quantised path integral, near the mass shell. But the regulator matrix M 
need not be diagonal. The off-diagonal terms of M encode the strength of an 'exponential 
of quadratic fiux'. These terms also depend on the choice of the time path, as we see from 
the difference in M for the contours C and Ci. 

The limit e ^ implies that the effect of the regulator matrix M is felt only on-shell 
(i.e. for — m? = 0). Equivalently, we may say that only world lines of infinite length 
(A = oo) see the regulator matrix. To see this, let M and M' be two different regulator 
matrices. We write 

L>M' =limi^oolim,^o+{[/ rfAe"*^^"^^^^ + / rfAe"*^^"^^^] 

Jo JL 

+ [ / c^Ae-*^^-^*^'^ - / rfAe-^^^-^^^] (3.5) 

J L J L 

iHX-eM'X / j\-iHX-eMX-i 



+ [ rfAe-^^^-^^^ - / c^Ae-^^^-^^^^]} 
Jo Jo 

The first square bracket on the RHS is Dm, the last vanishes with the indicated limits, 
while the second has support only on world lines of infinite length. Let us take M = I, 
the identity matrix. Then Dm = diag{ ^^2117^2^:^5 p2_^2_i^ }- Thus other propagators 
(obtained by using M not the identity matrix) differ from this basic one only by the 
contribution of world lines of infinite length. 

3.2. Quantising the relativistic particle. 

What is the origin of the two components a = 1 , 2 of the state on the world line? We 
would like to offer the following heuristic 'derivation' as a more physical description of the 
matrix structure in ( |3.2| ). 

The geometric action for a scalar particle is 

S = f mds = [ * m{X^',rX^,rf''^dT = [ Ldr (3.6) 

JXi JXi J 
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where r is an arbitrary parametrisation of the world hne. The canonical momenta 



satisfy the constraints 



P^P^-m^ = (3.^ 



We choose the range of the parameter r as [0, 1]. Following the approach in [|T^, we impose 
the constraint at each r through a 5-function: 



4n _ _ 

The path integral amplitude to propagate from Xi to Xf becomes 

G{X,^X,) = NJ y^ipig] e Jo (3.10) 

where is a normalisation constant, P^X^,t-= m(X^,T- )^/^ is the original action 
( ^■6| ) and we have divided by the volume of the symmetry group, which which is related 
to r-diffeomorphisms in the manner discussed below. (The 5-function constraint on the 
momenta and dividing by Vol[Diff] remove the two phase space co-ordinates redundant in 
the description of the particle path.) 

There are two ways to consider the symmetry of the action ( |3.1C1| ). The action is 
invariant under 

51 : SX^{t) = h{r)P^{r) 

5P^{r) = (3.11) 
5\{t) = h{r),r 

and 

52 : 5X^{t) = e(r)A(r)P^(r) 

5P^(t) = (3.12) 
d\{r) = (e(r)A(r)),. 

The difference between Si and ^2 is best seen by considering the finite transformations on 
A: 

Si : X'{t) = A(r) + ^ (3.13) 

52 : A'(r) = ^(r)A(r) (3.14) 
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Using Si we can gauge fix any function A(r) to any otfier function Ai(r), provided A, Ai 
have the same value of ^ 

A(r)dr = A (3.15) 



With such gauge fixing we obtain not a Green's function but a solution of the homogeneous 
Klein-Gordon equation. Thus the Fourier transform of ( |3.10| ) gives 

G{p) = 6{p^-m^) (3.16) 

for a suitable normalisation of the measure. 

With 1S2, A transforms as an einbein under the diffeomorphism r ^'i^)- Note 
that for regular e(r), A either changes sign for no r or for all r. We take Diff as the 
group of regular diffeomorphisms connected to the identity; then A does not change sign 
under the allowed diffeomorphisms. These diffeomorphisms cannot gauge-fix A(r) to any 
preassigned function Ai(r). We have the obvious restriction coming from ( |3.15| ), where A 
may now be interpreted as the length of the world line. This restriction is usually assumed 
to mean that the length of the world line is the only remaining parameter after gauge- 
fixing. But what we find instead is that there is a further complication: there is a discrete 
infinity of classes, each with one or more continuous parameters. One class comes from 
configurations A(r) which are everywhere positive. This class can be gauge-fixed with the 
allowed diffeomorphisms to have 

A(r) =0, / drXir) = A (3.17) 
Jo 

with < A < 00. Similarily, the set of everywhere negative A(r) can be gauge-fixed as in 
( p.l7|) but with — 00 < A < 0. Thus these classes give for the Fourier transform of ( |3.1C1|) 



Ga>o{p) =^j^ d\{r)e 



■i\/2(p'^(T)-m^) 



(3-18) 

Ga<o{p) = ^ [ cfA(r)e-^V2(p^(r)-^^) 



47r J_oo ^ ' p"^ — m? — ie 

respectively. (We have added the e term for convergence to each sector; this term was not 
in the action.) Keeping the first class alone gives the Feynman propagator for particles, 
while the second gives its complex conjugate.il Note that these two objects are the diagonal 
entries of ( ^.15| ) for /3 = 00. 



^ A restriction to the range (0, 00) for A can be naturally obtained using a Newton- Wigner 
formalism [ p^ . Here the particle travels only forwards in the time co-ordinate , thus it is not 
a CO- variant approach. 
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But in the path integral over A in ( p.lO|) we also have the class of A(r) which are 
positive for < r < ri, negative for ri < r < 1. The group of orientation preserving 
diffeomorphisms can gauge fix this to 



A(r) = for r ri, l drXir) = Ai, [ rfrA(r) = A2 (3.19) 

Jo Jti 

with < Ai < 00, —00 < A2 < 0. We would like to identify this sector as the contribution 
to the amplitude to start with a state of type 1 and end with a state of type 2 (the off- 
diagonal element D12 of the matrix propagator). A general sector has a given number 
of alternations in the sign of A, and in each interval of constant sign we gauge fix A to a 
constant A^. We can perform the integration over the variables A^ appearing in each sector, 
but we should also specify a 'transition amplitude' for each point Ti where A changes sign. 
We allow this amplitude to depend on the states on both sides of r^, and also on whether 
A changes from positive to negative or vice versa. This amlplitude is not supplied by the 
original action; it is supplementary information needed for determing the propagator. 

We would like to identify A > with particles of the type 1 field and A < with 
particles of the type two field in the language of section 2. To obtain the matrix propagator 
element Dji we need to add together all sectors for A(r) beginning as type i and ending 
as type j. Take the example (|3.2|), ( ^.3| ) and write 

/"OO _ /"OO 

D{p) = = / rfAe-'^^-^^^^ = / rfAe-^^^'-^^^^' (3.20) 
Jo Jo 



[p"^ - m'^ + ie{l + 2n{p))] 

- - ze(l + 2n(p))] 

-2^/n{p){n{p) + 1 



H' 

-2^n{p){n{p) + l) 
Expanding ( |3.20| ) in a perturbation series in M', 



(3.21) 



/"OO /"OO 

D{p) = / dXie-'^'^^ + / d\id\2e-'^'^^{-tM')e-'^'^^ + 
JO Jo 

= + D\-eM')D^ + ... 



(3.22) 



where 



-5 — z::^ — 7 



19 



Di2 for example gets a contribution from every term with an odd number of M' 
insertions: 

Di2(p) = -eD'.Mi^D", - £=fl?iM;2fl°2M^iD°iM;2D°2 + ... 

^ M»2 - + e2 ^ ^ 1 + 2n(p) ' ^ ' 

m=0 

= 2'k5{p^ — m^)\Jn{p){n{p) + 1) 

where e = e(l + 2n{j>)). 

Thus suppose the propagators for A > 0, A < are p2_^2-|-jg respectively. Let us 
choose the amplitude for A to change sign by requiring that momentum p goes to p, and 
the associated amplitude is '^^ ^ '^^i^^2n{p)^^'^ ~ ^sech(/3|po|/2). Then summing all the sectors 
that contribute to each Dji we get the matrix propagator ( |2.15| ). 

To obtain ( |2.12| ) , we must take the amplitude for A to go from positive to negative as 
esech(/?|]5o|/2)e^P°/2^ anti for negative to positive as esech(/3|po|/2)e~''^P"/2_ 

Note: To achieve the interpretation of a matrix of propagators as arising from the 
above 'sum over sectors', the matrix itself must be very special. From the symmetry of 
( p.l5|) we find that M must have the form ^ ^ ^ ^ . Because of the limit e ^ oo only the 

ratio h/a affects the obtained propagator. As we saw above, this ratio (which directly gave 
the amplitude for reversing world line orientation) got determined by the value of D12. 
But now there are no free parameters to adjust for obtaining Dn. A direct calculation 
gives 



00 ,Vn~V^MMPl±11^2m-n2m 



^ ^(^^ l+2n(p) 



{P^ '~ + ^^)((p^ ~ m?)"^ + e^)"^ 
^ + 2ti5{p^ -m^){l + 2n{p)) 



(3.25) 



p2 _ jy^2 _|_ 



which agrees with Dn in ( |2.15| ). Thus the matrix of propagators obtained from the real 
time formalism has, at least for this example, the special form required for the first quan- 
tised interpretation to work. ( ^.121 ) also has the required form. M has unequal off diagonal 
terms but their product remains the same as in the previous example; thus we get the same 
Dii, in agreement with ( p.l2|) . 

Thus we have verified in these examples that the amplitude for einbein sign change 
reflects the presence of particle flux in the first quantised formalism. In section 4 we will 
present an example where the flux created by spacetime expansion is handled in a similar 
way. 
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3.3. BEST formalism. 

One might wonder if the more formal BRST quantisation of the relativistic particle 
would resolve the above issues about different possible quantisations. We follow the nota- 
tion in 0]. We introduce the canonical conjugate tt for A ([A, tt] = i) and ghosts (rj^^Vi), 
{'rf:V2) = foi^ the two constraints tt = and ^(-P^ — m^) = respectively. 

The BRST charge 

2 

g = r^^TT + —{P^-m'^) (3.26) 

is nilpotent (Q^ = 0), and gives the variations: 

SX^ = ~i[X^, Q] = rfP^ 5Pf, = 
SX = r]^ 6n = 

5r]^ = 5Vi = -Tx 



(3.27) 



The equation of motion (obtained after gauge fixing) gives t)^ = i]^, which agrees with 

But we can define another nilpotent BRST charge 

Q' = r]^'n'X' + %'2(p2_^2^ (328) 

which generates the symmetry 

SX^" = r]^'X'P^' SP^,' = 

6X' = X't]^' 6n' = —T]^ tt' ~ 

dr]^' = 5V[ = -X'n' ~ 

V = SV^ = --{P^' -m^)^0 

2 



(3.29) 



The equation of motion gives rf' = r]^' X\ which suggests that we should identify the above 
symmetry with 82- 

The symmetries Q and Q' are related through the identifications 



TT 



tt'A', a = logX' (3.30) 



all other primed variables equalling the unprimed ones. From ( p.30|) we find that — 00 < 
A < 00 corresponds to < A' < 00. If we perform a path integral with the primed variables 
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and sum over both positive and negative A' then we are summing over more than is being 
summed in the unprimed variable path integral. 

We thus see sources of ambiguity on the quantisation of the relativistic particle work- 
ing with a Fadeev-Popov approach in sec 3.2 and a BRST approach in sec. 3.3. In fact 
the action we start with, ( p^ ) , is itself ambiguous because of the two possible signs of the 
square root. The particle trajectory would keep switching in general between timelike and 
spacelike, and at each switch we have to choose afresh the sign of the real or imaginary 
quantity obtained in these two cases respectively. This suggests that the world line con- 
figuration should be described by the pair {X^{r), cr(r)} with a = ±1 giving the choice of 
root. Evaluating the quadratic form of the action (given in (|3.10|) ) classically we find the 
sign of A to be related to the sign of the square root chosen for (^.61). 

The above discussion suggests a close connection between the ambiguities found in 
three different approaches to the quantum relativistic particle, and it would be good to 
determine if they indeed are the same. For the rest of this paper we simply adopt as 
basic the picture of two complex conjugate propagations on the world line, with switching 
between them possible through the regulator matrix. 

3.4- Summary. 

We may describe our quantisation of the relativistic scalar particle in the following 
colloquial terms. To obtain the usual Feynman propagator the particle was allowed to 
move both backwards and forwards in target space time. But it moved only forwards in 
its proper time. Now we are allowing the particle to move both backwards and forwards 
in proper time also; this corresponds to flipping between the two signs of the world line 
einbein.i Giving the amplitude to reverse orientation of proper time amounts to specifying 
an 'exponential of quadratic' density matrix,. Putting this amplitude to zero lets us recover 
the Feynman propagator, by restricting to the positive sign of the einbein. 

Thus we find a covariant description of the notion of the notion of density matrix. 
We see that in the first quantised dsecription, there is very little difference between causal 
perturbation theory for pure states and for mixed states. 

Note: trp is finite only for Aaj < (1 — e~^)^ (see appendix). Consider the initial 
value problem for charged scalar field in a small box with a (suddenly switched on) strong 
electric field. We may wish to expand the field in modes which are eigenfunctions of the 

^ I am indebted to A. Vilenkin for this way of expressing the result. 
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Hamiltonian. The above inequality is violated at the point where the electric field becomes 
strong enough to destabilise the vacuum PO] . For stronger fields the mode coefficients are 
6, d^, [b,d^ = — z instead of a, a\ [a, a^] = 1- We expect that black holes would also be 
characterised by a similar violation. The density matrix is still 'exponential of quadratic 
in the field' but (|2.25 ) is not a useful representation. We will ignore this limitation of the 
form (|2.25| ) because for our examples it will be adequate. 



4. A curved space example: spacetime with expansion. 



Consider the free scalar field (( p.l| ) with A = 0) propagating in 1 + 1 spacetime with 
metric 

ds"^ = C{r])[dr]'^ - dx'^], -oo <?7<oo, 0<x<27r (4.1) 

C{r]) = A + Btanhnr], A>B>0 (4.2) 

The conformal factor C{r]) tends to A ± B at r] —>■ ±oo. The limit k ^ oo gives a step 

function for C{r]); the Universe jumps from scale factor A — B to A + B at r] — 0. We will 
work in this limit to ensure simpler expressions. 
For ?7 ^ — oo it is natural to expand cj) as 

°° 1 1 

(^(?7,x) = y ^^=(a^e-*'""''+*"" + aj.e^'"" "-*"") (4.3) 

with 

uj- = (n2 + (A-S)m2)i/2 >o (4.4) 
We define the 'in' vacuum by 

flnlO >in= 0, for all n (4-5) 

Similarily, for r] — > oo we write 

oo 

with 

0J+ = (n2 + (A + S)m2)i/2 >0 (4.7) 



23 



The 'out' vacuum is defined through 



an\0 >out= 0, for all n {4.t 



The 'out' vacuum does not equal the 'in' vacuum, even in the free theory: 

n>0 



|0 >out = Coe'-^''o4 Yl C^e'^-^^^-HO (4.9) 



C„ = (1-62)1/2, ^>0. (4.10) 



As mentioned in the introduction, evaluating the analogue of the Polyakov path inte- 
gral gives (we denote the pair (77, x) by z) 

JO out < > 

(We have integrated out p in the phase space path integral.) It is not surprising that both 
vacuua appear in this quantity; after all the action and measure are covariantly given and 
should not distinguish the past or future as special. Using —i in place of i in the exponential 
gives in < 0\f[(p{z2)(p{zi)]\0 >out I in < 0|0 >out- {f deuotes anti-time-ordering.) 

We wish to compute the propagator appropriate to the initial value problem based on 
the 'in' vacuum at t = —00. This implies a density matrix p = po = |0 >m in < 0|. To 
develop a perturbation theory using po we need a real time contour running from 77 = — 00 
to 77 = 00, and then back to 77 = —00 where we insert po and take a trace to close the path. 
This perturbation theory requires a matrix propagator D: 

r,f, ,N _ f ^n < 0\T[(P{z2)(p{Zi)]\0 ^n < 0\(p{Zi)(P{z2)\0 >^n \ . ^^ 

^ V ^n<0|(/>(^2)</>(^l)|0>.. <0|T[(/>(^2)</>(^l)]|0>.J ^ ^ 

Our goal is to see if there exists a choice of regulator matrix M such that evaluting 

^ - X ' -iU + m^)) (4-13) 

gives the matrix propagator ( [1.12| ). 

Let us discuss the general structure of the world line Hamiltonian and regulator matrix 
in ( ^.13| ). H and M have a '2x2 matrix of operators' structure, on account of the two signs 
of the einbein A possible on the world line. The operator □ + m? acts on the space of 
all functions of (77, x). Formally, each operator in the 2x2 matrix M also acts on this 

24 



space, but as noted at the end of section 3.1 the only space affected by these operators 
is the space of 'on sheU' wavef unctions, i.e., solutions of (□ + m?)f{ri,x) = 0. In the flat 
space examples ( |3.ij| ), ( ^.4| ) we had translational invariance in space and time, so each 
operator could connect one Fourier mode p only to itself. Thus for each p, M was just 
a 2x2 matrix of numbers. In the present example, we retain x-translational invariance, 
so different x-Fourier modes (labelled by n) are not mixed by M. But time translational 
invariance is broken by the expansion, and so the two solutions of the wave equation for 
fixed n can be transformed into each other by the elements of M. Thus for each n we 
need to consider a 4x4 matrix, which acts on a column vector {{fn: fn}~^ ^ {fn^ fn}~)- Here 
the first pair of functions propagate on the world line as e*^'-'"'""^ ^ while the second pair 
propagates as e"*'^'-'"'""^ \ Within each type (+ or — ) we have two linearly independent 
solutions of (□ + m?')f = 0. Thus for example the element M41 of the regulator matrix 
gives the amplitude for the einbein to change sign from positive to negative, and for the 
wavefunction on the world line to change from to 

Let us set up the calculation of Green's functions in the first quantised formalism. We 
need eigenfunctions of the world line Hamiltonian: 

(□ + m^)?/^s(?7, x) = -Hi/jsiv^x) = si)s{r},x) (4.14) 
The following is a complete set: (— 00 < n < 00) 
2 , ^ ^ 2 



m H — — ■ < s < m + 



A+B A-B 



27r 

^"^[1(1 + Z+ )e— + 1(1 _ ^)e''^-^], ry < 



V27r 2 w_ 2 iu- 

i>+ = [-{A + B){m^ - s) - n'^]^/'^ > 0, u_ = [{A - B){m^ - s) + n^]^/'^ > (4.15) 



-00 < s < + 



n 



2 



A + B 



V ZTT 



^"^[1(1 + ^)e— + 1(1 - ^)e^"-^], 77 < 



/27r 2' 2' u. 
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These functions are normalised as 



sign(z^_l_) = sign(z^_) 



(4.16) 



2 ~2 



(4.17) 



/ \ 2 2 2 

(V'.' n'^V'.+.n) = 5r.>,nAA + B)[^^^-^^5iu'+ - u+) + + u^)] (4.18) 

Let us first recover the propagator (|4.11j ) in this formahsm. Let ry', ry > 0. The range 

gives the contribution 



py/2B/(A-B)\n\ roo 







<ilJi>'^,n\v.x> d{u+-i>+)[{A + B)Tv] ^^^^-^ 



v_ + v_ 



+ 



1 



(4.19) 



Similarily the range —00 < s < w? + provides the contribution 



E 



i(x'— x) /-oo 



27r 



dv 



H ^ e" 



1 



^ {yl - n2 - w?{A + B)+ it) 

(4.20) 

There is a branch cut in the complex z/_|_ plane joining z/_|_ = ±«i/ ?_^d |n|. The vj^ integral 



in (|4.20|) has a discontinuous jump across this cut for the part multiplying e^^+^'^ The 
contribution from p.l9| ) can be added to ( |4.20| ), however, with the identification z?_|_ = 
This results in a contour passing over the cut. Evaluating the resulting contour integrals 
one obtains the result 



^in{x' —x) 

2ti 2ujt ' 



2UJ n 



i^+iv'+v)^ for 77,77' > 



(4.21) 

which may be readily verified in the operator language using (|4.6|) and ([4.9| ). 

We now proceed to the 'in-in' matrix propagator. From the expansion of the field 
operator (p in creation and annihilation modes, we can readily compute each element of 
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the matrix propagator. We try to find a regulator matrix M which wiU reproduce this 
propagator by the first quantised calculation. To display the 4x4 matrix M, we need to 
choose a basis in the space of solutions of the wave equation for each x-Fourier mode n. 
We choose the basis {u~ > 0) 

= ^e-^-" ^ rj < 



'2n 



v^27r '■ 2 ^ UJn 2 u;n 

The result we get is that the required propagator is obtained for 

/I q \ 

M = 



1 

-2 4:Br 

^4B^ 1+41 





where 











-2 




1+4B 


1 








1 







(4.23) 



Bn = -^^^-^ (4.24) 

^ UJn + UJn 

(In this computation we need to note that there are many equivalent descriptions of the 
delta function. For example lime^o ^2^^2 = 7r5(x), but limg^o (x+ie^ ~ ^' ^^^^ amount 
of the latter can be added to the former to get a different representation of the 5-function.) 

As was the case in the examples of the last section, M does not have enough indepen- 
dent entries to reproduce any arbitrary 2x2 matrix of propagators. In the present example 
any propagator has the freedom of 4 constants for each x-Fourier mode n: we can add 
Ylf j^i Cij f^^\z) f^^\z') where Z*^*-* are given in ( ^4.22|) . A 2x2 matrix of propagators thus 
has 16 free parameters for each n. But the entries in the diagonal 2x2 blocks of M are not 
all independent, in the sense that some changes to M do not affect the final propagator 
obtained. (For example, scaling M cM leaves the propagator unchanged.) Thus at first 
it seems a coincidence that we could find M to obtain (|T|). 

To understand this coincidence we recall that ( [4.12|) came from a 'closed time path' 
formulation. This formulation is really covariant (and not dependent on any choice of 
'time'). Each spacetime point is covered twice in the description (alternatively we can say 
that there is a doubling of fields ) . Over this double cover of spacetime we are looking 
for a first quantised description of the propagator. Clearly such a propagator can be 
represented by the usual path integral in each cover, together with an amplitude to switch 
between covers. Thus all matrix propagators arising in the closed time path formalism 
will have the special form that arises from allowing orientation changes on the world line 
in the first quantised language. 
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5. Strings. 

Let us now come to strings. Requiring consistent propagation of the string determines 
the field equations for the background on which the string moves. But use of 'in-out' and 
'in-in' propagators might give diffrerent field equations, as we argue in this section. 

To obtain the necessary ingredients for the argument, we first consider a simple case: 
strings in a contant temperture bath, in a time independent background. Since we are 
at constant temperature, we can use for the moment the periodic imaginary time trick to 
take into account the temperature p^ . 

We expect that the vanishing /3-function conditions will reproduce, in some approxi- 
mation, the Einstein equation 

C^. =<T^u> (5.1) 

where T^j^ is the stress-tensor of the thermal distribution. The tree level /3-function cal- 
culation is a local one: local on the world sheet and local in target space. At tree level 
we thus obtain G^^, = 0; we are unable to 'see' the thermal bath. The /3-function also has 
a contribution at one loop in the string coupling expansion: we must consider tadpoles 
attached to the world sheet by necks thinner than the world sheet cut-off scale [|3| . The 
loop of the tadpole can wind around the compact imaginary time direction. This takes 
into account the value of temperature and produces the contribution < T^,^ > required for 
the field equation |2|]. 

Let us now ask how the /3-function equations would reproduce (|5.1| ) if we are not 
permitted use of the periodic imaginary time trick to take into account the particle distri- 
bution. The question is pertinent, because as mentioned in the introduction time indepen- 
dent distributions are unnatural in a theory of gravity, and any time dependence would 
invalidate the periodic imaginary time trick. There is no 'analytic continuation of time' 
for a general geometry. @ 

In the absence of the imaginary time winding mode, neither the tree nor the loop 
contributions to the /3-function see the thermal distribution. But considering the extended 
propagator-vertex structure of the real time formalism we should be able to recover the 
RHS of The one loop tadpole is composed of two propagators and two vertices 

(including the vertex joining the tadpole to the rest of the world sheet). These vertices 



^ It is possible to study some time dependent situations by making a canonical transformation 
on the target space co-ordinates and taking the new X° to be imaginary and compact; this 
does not however have the physical appeal of the real-time method. 
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can be type 1 or type 2, in the notation of section 2, and the propagators wiU be the 
corresponding elements of the 2x2 matrix of propagators. These propagators depend on 
the ensemble distribution function, and so 'see' the temperature. It would be interesting to 
explicitly carry out this calculation, but here we will just assume that the real time method 
will work for strings the same way it does with particles. Note that such a calculation does 
not need the restriction to time-independent distribution functions. 

We now argue, again with the scalar field analogy, that there is a direct connection 
between the choice of propagator and the value of < T^^ > seen by the /3-function calcula- 
tion for the string background. Consider the scalar field in flat space and ignore the mass 
term for the moment (since the string case has massless fields on the world sheet). Then 
Einstein's equations read 

= < d^<P{x)d,<Pix) > = lim d^dl < (P{x)<Pix') > (5.2) 

X — >x' 

(We ignore any subtraction in the definition of T^i, since the string case is ultraviolet finite 
and has no subtraction.) Using the analogue of the Polyakov path integral 

(5.3) 

where ds = Xdr is the length element along the world line. Then 

a,al < 0(x)*,(.') > = ^"dAflmJJJ);:' ^(0)^^(A)e-^r """^'-'S (5.4) 
Let us compare the above with the /3-function equation for string theory in a finite tem- 



perature bath described in the imaginary-time formalism 

2Til3Vd 



2cj> 

= i?^. + 2V^ V, = ^^^j^ Yl < 9X^dX^ >g (5.5) 



[(j) is the dilaton, V is the voloume of spacetime, d is the number of space dimensions.) 
The RHS of ( |5.5[ ) is the analogue of (|5.4|). (The factor {I3V)~^ in ( |^ ) compensates for 
the X^ zero mode; without this compensation we would compute the stress tensor not at 
a point X but instead integrated over all the range of coordinates X^.) Examining the 
derivation of (|5.4[) , we expect that in the string case also the stress tensor contribution to 



the /3-function comes in a direct way from the propagator of the first quantised language. 
(Thus in ( ^.5| ) it comes from the propagator in spacetime with a compact Euclidean time 
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direction.) If we have a general time dependent spacetime (with Minkowski signature) and 

we compute amphtudes by the Polyakov prescription, then the propagator is 

out < 0|T{^(x)(/)(x')}|0 >in I out < 0|0 >in and so we we expect to incorporate 

out < 0\T^u\0 >in I out < 0|0 >in in the background field equations obtained from the 

/3-function constraint. This is the wrong kind of object to have in the gravity equation, as 

we illustrate with the following simple example. 



We consider the free scalar field on the 1 + 1 spacetime (|4.1| ), ([4.2|) . We compute 'in- 
out' vacuum expectation values and 'in-in' vacuum expectation values, for the stress tensor. 
To achieve a parallel with the string case, we consider no subtraction in defining T^^, but 
compute the contribution from different field modes separately, to get finite answers. (We 
can imagine an ultraviolet cutoff regulates the sum over different modes, to make sense of 
sum over modes.) More specifically, we compute the contribution to 

= / dx in < OlT^iy^T], x)\0 >iri 
Jx=0 



r.27r 

in out = I dx out < ^\T^i/{r],x)\0 >in /out < 0|0 >i 

'x=0 



(5.6) 



from the field modes with a;-Fourier component ±n: 

a±n f±n + a±n f±n (5-7) 



(P±n 
p±inx 

_ e^^"'" ^ u++uj- ^ UJ+ - uj- ^^^^^^ ^^^^ 



(5.8) 



The stress tensor for the scalar field is 

T^. = lidf,cf)d,(j) + d,(j)d^(j))-^g^Md^(p+^m''(j)'' (5.9) 



T^Af^a] = ^{djd,g + djd^g)-^g^,dxfd^9 + ^m^fg (5.10) 



Define 



The expansion ( |5.7| ) is in terms of the 'in' vacuum creation and annihilation operators. 
We recall that 

|0>o.t = Ce'"""-"|0>^n (5.11) 
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where b = 



— (cu-i- — a;_)/(a;_|_ + u;_). We have 



< Too >in in = Too[fn, /*] + Too[/-n, f-n] 

< Til >in in = Tii[fn, f*] + Tii[/_^, /*^] 

< Too >in out = Tooifn, fn] + ^Oo[/-n, f-n] + ^{^Ool/ni fn] + ^Oo[/-n5 f-n]} 

< Til >^n out = Tii[fn, f:] + Tn[/-n, fin] + b{Tii[f:. f:] + Tii[fl^, 

(Toi, Tio vanish.) 

From (^]^) we get 

< Too >in in = V <0 

< Til >in in = n^/uj-, rj <0 

2 2 

/ 2 2 \ / 2 2 \ \ r\ 

= - ^{uj^+uj_) - —^{u^-u_)cos{2u+t), r)>0 

ZUJ—UJ_^ ZUJ_^ 
< Too >in out = UJ-, r] < 

= 1] > 



(5.12) 



(5.13) 



71^ (a;_ - n^) uj+ 



< Til >^n out = —- r:^_Z(r:^_i^)e2--'', ^ < o 

00- 00- OOj^ + U- 

= Z^_(^4^(^^^±^)e--.., ,>0 

Here < Nout >= ('^+ — u;_)^/4c<;_(_c<;_ is the number of 'out' particles in the 'in' vacuum. 

We see that < Tqo >in in has a direct physical interpretation. For ry < (before 
expansion) we are in the 'in' vacuum state, and we just get the vacuum energy oo-/2 for 
each of the two modes n, —n considered here. For 77 > (after expansion) we get the 
vacuum energy ooj^/2 for each mode, plus the energy of created 'out' particles. 

< Too >in out on the other hand gives only the vacuum energies; there is no contri- 
bution from created particles. Worse, < Tn >in out is complex, and so cannot be on the 
RHS of Einstein's equation for the gravitational field. 

This result is not surprising; we know that in the 'in-out' formalism the effective field 



equation for a real field can become complex [|Tl|. The 'in-in' formalism, by contrast, 
is causal and yields real expectation values for self-adjoint operators. Thus it appears 
essential to extend the first quantised path integral prescription for strings so that we can 
get the 2x2 matrix propagators of the real time formalism. Extending the discussion of 
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section 3 to the string case we would obtain a sum over both signs of the zwiebein on 
the world sheet (instead of just a sum over metrics.) The amplitude to flip sign of the 
zwiebein would encode a density matrix which describes 'exponential of quadratic' string 
distributions (initially present or 'created'.) Because the string is two dimensional, there 
is an interesting variety of ways that the zweibein can change sign, but for the purposes 
of the present paper we can just assume that the string is expanded into particle modes, 
in which case the above approach for particles applies. 

Note: Leblanc I^S]] studied the real time formalism for open and closed strings, for the 
case of constant temperature in flat space. The propagator was computed in the 'thermo- 
field dynamics' language, which used the Niemi-Semenoff time path, and so was given 
by (|3.2|) , (|3.3|) . For this time-independent situation amplitudes were computed and the 
Hagedorn temperature recovered. To extend the string calculations to arbitrary spacetime 
geometry, as required in our approach, needs further progress in defining and computing 
amplitudes in Minkowski signature spacetimes. We will discuss this issue in more detail 
elsewhere. 



6. Discussion. 

Let us summarise the arguments and results of this paper. Requiring consistent prop- 
agation for the first quantised string should determine the background geometry. Do the 
background field equations incorporate the stress tensor of 'created strings' in time depen- 
dent geometries? To investigate this question we studied a simpler theory in first quantised 
language: the relativistic scalar particle. From the example of strings at constant tempera- 
ture, we expect that backreaction from particle fiuxes will show up at one string loop level 
in the /3-function calculation. But this one loop contribution computes the stress tensor 
as two derivatives of the propagator, and so the backreaction seen by the /3-function cal- 
culation depends on the choice of propagator. The Polyakov prescription of summing over 
target space co-ordinates and world sheet metrics gives an 'in-out' vacuum propagator. 
Using the scalar field analogy, we then argue that the gravity field equations obtained at 
one loop level will have an object like out < 0\Tfj,u\0 >in /out < 0|0 >in for the backreac- 
tion. We demonstrated by a simple example that such quantities are in general not real, 
and not the correct ones to have in the field equation. We need instead a 'true expectation 
value' of the stress tensor, which will be achieved if we have a causal 'in-in' propagator, 
instead of the 'in-out' propagator. 
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In achieving a causal formulation of perturbation theory we need to extend the prop- 
agator to a 2x2 matrix of propagators [|Tl[] . This closed time path formulation (called the 
real time formulation in non-equilibrium many body theory) also naturally handles per- 
turbation theory in the presence of particle fluxes. In fact the perturbation scheme is built 
around precisely the same kind of fluxes (exponential of quadratic) which are obtained in 
particle creation effects. 

Accepting the necessity of a real time formalism for strings, we are faced with the 
question: can this formalism be obtained in a natural way in a flrst quantised language? 
We start with the square root action of the scalar particle, and try to reach the quadratic 
form that is the analogue of the Polyakov action. We flnd that the proper time in the 
latter description does not necessarily go forward. We need to supply the amplitudes for 
any state > on the world line to connect to any other state >, when the proper 
time reverses direction. Setting all these amplitudes to zero is a special choice that forbids 
reversal of proper time; choosing the forward direction of this time reproduces the 'in-out' 
vacuum propagator. 

But it is also possible to set these amplitudes to be nonzero. With a particular choice 
of the amplitudes we found the matrix propagator for perturbation theory in a thermal 
bath. With another choice we found the causal 'in-in' propagator for a Universe which is in 
the 'in' vacuum but where expansion creates a bath of 'out' particles. In both these cases 
the real time description involved a density matrix of the kind ( |2.25| ). For the thermal 
bath we have a = 7 = 0, /9 flnite. For the 'created particles' case we had a limit for the 
coefficients in which a pure state density matrix was obtained. 

For strings, the passage from the Nambu-Goto action to a quadratic action would 
bring in two signs of world sheet zweibien; restricting to one sign would give the Polyakov 
amplitude. As remarked at the end of section 3.1, only infinitely long world lines contribute 
to the orientation fiip amplitude of the world line. The analogue of this for the string case 
would be that the effect of p is felt only through the boundary of the moduli space of 
Riemann surfaces, where a homologically trivial or non-trivial cycle is pinched. A f3- 
function calculation for the string world sheet theory would have to take into account such 
pinches while considering the small handle contribution studied by Fishier and Susskind. 
Such a calculation would yield a relation between the classical fields and the particle 
fiuxes, rather than just among the classical fields giving the background. (The constant 
temperature case quoted in section 5 is an example of this.) The classical field strength 
describes a coherent state; the fiux describes an 'exponential of quadratic' distribution of 
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particles. The natural way in which field and fluxes have been mixed in our formalism 
leads us to expect that in a theory of gravity 'effective equations' of the quantum theory 
should have both these components, instead of having field values alone. 

It would be interesting to obtain a non-perturbative treatment (like a matrix model 
approach) for a situation like an evaporating black hole where the backreaction from a flux 
is important in determining the evolution of the geometry. Based on the discussion of this 
paper, it appears that for this problem we should be looking not for a conformal theory 
but an 'extended conformal theory' (with both orientations of the world sheet allowed). 
The conformal invariance condition would set up a relation between the coherent field 
(contributing to the /3-function from tree level onwards) and the orientation reversing 
amplitude representing 'string fiux' (which contributes to the function from one string 
loop onwards). The conformal theory describes just one string, and so should be an 
inadequate description of this situation. 
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Appendix A. Wick theorem for 'exponential of quadratic' density matrices. 

We wish to establish Wick's theorem for density matrices of the form 

A string of creation and annihilation operators can be brought to normal ordered form 
in the same way as for the usual Wick theorem in the vacuum. What we need to show in 
addition is that 

J-TV{pat...ata...a} = Y ^^^^'^'^ . . . ^^^^'^^ . . . ^^^^^^^ (A.2) 

where the RHS has a summation over all possible pairings of the a^, a operators on the 
LHS. We sketch below some of the steps involved in the derivation. 
Note 

^aa'^ a'^ ^-I3a^ a _ ^- /3a^ a ^a' (A3) 

with a' = ae^^. A straightforward calculation gives 



4a7 

which we may rewrite as 



= (l-^"')"^[l-(Tf^]"^/^ (A.4) 



Trp = Tr{e^""e-^"^"e"'"^"^ = (l-e-^)-^!- , ^""X^J "^^^ (A.5) 

(1 — e P)^ 

From ( |A.1|) we see that in (|A.2| ) there must be either an even number (2p) of 'a' oscillators 



and an even number (2^) of 'a^' oscillators, or an odd number (2p + 1) of 'a' and an odd 
number (2g + 1) of 'a^' oscillators. Assume first that we have the former case. Then the 
LHS of (|A.2[) is obtained as 

Tr{p(at)(2'?)(a)2?'} = A^^d^,)<i^d^)PTrp (A.6) 

Trp 

(Partial derivatives are taken with a', /3, 7 as independent variables, unless otherwise 
mentioned.) In particular, 

< a^a^ > = ^a^,Trp = 2-fe-'^^/K = A 
Trp 

<aa> = T^djTvp = 2a e'^^/i^ = B (A.7) 
<a+a>= - ^c>^[IVp]„,^ = e-^{l-e-^)/K = C 
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where in computing C, dp is a partial derivative with a, 7 held fixed. Here 



K = (l-e-^)2-4a'7e"'^ (A.8) 

We find 

Trp = e^[C2-AS]i/2 = K'^''^ (A.9) 

For fixed /3 



do,' A 


= 2A^ 


d^A = 


2C^ 


do,'B 


= 2C^ 


d^B = 


2B^ 


do^'C 


= 2AC 


d^C = 


2BC 



(A.IO) 

Using the above formulae, we can establish (|A.2|) by induction. Suppose (|A.2|) holds 



with 2p operators 'a' and 2q operators 'a^'. A typical term on the RHS would have the 
form FA^^B^^C^^, where F is a constant and ni, n2, ns > 0. To establish the result for 
2p operators 'a' and 2q + 2 operators 'a^' we get for the LHS of ( |A.2| ): 

d' [e^fC^ - ABV^'^FA'^^ S"2C"^3l ^ ^Mni + l^na^ng 

+ 2?iiyl"i+iS"2C"^3 + 2?i2A"^S"2-iC"3 ^ 2n3A'^i+^S"2C""3] 



The first term on the RHS of (|A.ll ) gives the pairing of the two new operators with each 



other. The second term gives the ni ways to choose an existing pair (a^a^) and to contract 
the new operators with members of this pair instead. The third term corresponds to 
choosing an (aa) pair in the original expression and contracting the 'a' operators with the 
new 'a^' operators instead. The last term corresponds to exchanging the in an existing 
a^a pair with one of the new operators. It is easily seen that this generates all the new 
terms required on the RHS of ( [A.2| ) for the induction to hold. 

To work with the case of an odd number of a and operators we start with the 
expression CTrp = Tr{e'''""e~'^"^"e" "^"V^a} in place of Trp, and proceed as above to 
introduce extra a^a^ and aa pairs in the induction. 
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